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FOREWORD 


This  Final  Technical  Rfport  covers  the  work  performed  under  Contract 
No.  F33615-72-C-2061,  Project  No.  7633,  from  15  June  to  15  September  1972. 

The  contract  with  the  Sperry  Rand  Research  Center,  100  North  Road, 
Sudbury,  Massachusetts,  01776,  is  to  investigate  several  techniques  for  im¬ 
proving  the  accuracy  of  an  Air  Force  Avionics  Laboratory  time  domain  measure¬ 
ment  system  for  measuring  the  constitutive  properties  of  radar  absorbing 
materials. 
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SECTION  1 
INTRODUCTION 


1.1  BACKGROUND 

The  investigations  described  in  this  report  were  performed  at  the 
Sperry  Rand  Research  Center  (SRRC)  in  Sudbury,  Mass,  for  the  Air  Force  Avionics 
Laboratory  (AFAL),  Wright-Patterson  Air  Force  Base,  Ohio  under  Contract 
F33615-72-C-2061  over  a  period  from  June  through  September  1972. 

In  September  1971  SRRC  delivered  to  AFAL  a  unique  prototype  system 
for  measuring  the  constitutive  parameters  of  materials  over  the  frequency 
range  0.1  to  10  GHz  without  disturbing  the  sample.  This  system  was  based  on 
time  domain  metrology  principles  where  the  transient  responses  of  a  material 
to  a  smoothed  impulse  serves  to  characterize  the  test  material.  The  basic 
components  are  a  subnanosecond  baseband  pulse  generator,  coaxial  delay  line 
networks,  a  sampling  oscilloscope,  a  controller,  magnetic  tape  recorder,  and 
analog-digital  interfacing.  The  time  domain  responses  are  measured  and  record¬ 
ed  on  magnetic  tape,  and  a  subsequent  Fourier  transform  program  yields  the 
desired  complex  permittivity,  e  and  complex  permeability,  n  ,  More  detail¬ 
ed  information  on  this  system  and  its  operation  is  contained  in  two  recent 

1  2 

AFAL  Technical  Reports. 

During  the  first  year  of  operation  at  AFAL,  the  system  has  proved 
to  be  extremely  useful  for  measuring  the  constitutive  parameters  of  lossy 
materials.  As  a  result  of  usage,  certain  accuracy  limitations  were  noted  in 
the  prototype  unit.  The  objective  of  this  program  was  to  investigate  several 
techniques  for  improving  the  accuracy  and  extending  the  capability  of  the 
AFAL  system. 

1.2  TASKS 

The  tasks  in  this  investigation  fell  into  four  categories: 

a)  The  primary  effort  was  a  coaxial  line  investigation  to 
reduce  the  observed  systematic  error  in  the  calculated 
values  of  e  and  p,  .  This  included  the  design  and 
fabrication  of  necessary  modifications. 

-  1  - 


b)  A  study  to  establish  the  cause  of  and  to  eliminate 
sporadic  low  frequency  errors. 

c)  An  investigation  of  the  feasibility  of  using  standard 
samples,  with  known  e  and  p,,  to  calibrate  the  system  and 
to  develop  software  to  impleme;.-  this  technique. 

d)  A  study  to  extend  the  capability  of  the  system  to  the 
measurement  of  low  loss  dielectric  materials. 

Each  of  these  tasks  is  described  in  subsequent  sections.  In  addition,  a 
preliminary  safety  analysis  was  oerformed  to  minimize  catastrophic  failure. 

Dr.  Harry  M.  Cronson  was  the  principal  investigator  in  this  program. 
Mr.  Peter  G.  Mitchell  was  responsible  for  hardware  modifications  and  perform¬ 
ing  experimental  measurements.  Technical  assistance  was  provided  by  Mr,  Richard 
L.  Earle. 


SECTION  2 

COAXIAL  LINE  INVESTIGATION 


2.1  INTRODUCTION 

The  major  effort  in  this  contract  was  the  reduction  of  systematic 
errors.  These  errors  have  the  appearance  of  ripples  on  a  plot  of  the  real 
and  imaginary  parts  of  the  relative  permittivity,  e  =  e'-je",  and  the  rela¬ 
tive  permeability,  p'-jp",  vs  frequency.  Unlike  ripples  caused  by  random 
noise,  the  peaks  and  valleys  of  systematic  ripples  usually  occur  in  the  same 
frequency  range  regardless  of  the  material  tested.  The  systematic  ripples 
are  usually  larger  than  random  ripples  and  cannot  be  reduced  by  averaging. 

The  probable  cause  of  these  systematic  rippies  is  mismatch  in  the  connector 
and/or  RG331  foasu  cable  in  the  vicinity  of  the  sampler  holder.  (The  AFAL 
10  ns  measurement  system  was  constructed  from  semi  flexible  RG331  cable  in 
order  to  achieve  a  relatively  compact  configuration.)  On  tne  other  hand, 
systematic  errors  were  not  observed  in  the  prototype  SRRC  2.5  ns  system  con¬ 
structed  from  14  mm  precision  air  line.  Since  the  connectors  to  the  foam 
cable  present  a  larger  discontinuity  than  General  Radio  (GR)  air  line 
connectors,  this  discontinuity  could  cause  the  observed  ripples. 

Our  method  of  approach  was  to  divide  the  coaxial  line  investigation 
into  4  sequential  subtasks: 

a)  Obtaining  a  quantitative  measure  of  the  magnitude  of  the 
ripples  of  resting  systems  for  comparison  with  later  results. 

b)  Determining  the  cause  of  the  ripples. 

c)  Investigating  methods  of  ripple  reduction. 

d)  Fabrication  3nd  testing  of  the  recommended  modification. 

The  remainder  of  this  section  describes  work  done  under  each  subtask. 

2.2  MEASURE  OF  RIPPLE  MAGNITUDE 

To  reduce  the  observed  spurious  responses  in  p  ano  5  ,  it 

seemed  reasonable  to  first  establish  the  maf-itude  of  these  ripples.  These 
numbers  could  then  serve  as  a  basis  of  comparison  for  an  improved  system.  As 
a  measure  of  the  ripples,  we  sought  an  easily  computed  quantity  related  to  the 
difference  between  the  maximum  and  minimum  observed  value  of  ,  The 
quantity  arbitrarily  selected  and  called  the  percentage  ripple,  R  ,  is 
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(1) 
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It  should  be  noted  that  the  observed  ripples  are  due  to  both  random 
and  systematic  errors.  However,  since  the  systematic  errors  are  usually 
larger,  the  above  equation  furnishes  a  useful  guideline  for  systematic  errors. 
To  avoid  extreme  points  at  both  ends  of  the  spectrum,  where  random  noise 
errors  are  large,  the  0.1,  0.2,  0.3  and  9.8,  9.9,  and  10  GHz  values  of  e' 

it 

were  not  included.  From  available  data  it  was  convenient  to  choose  a  Teflon 
sample  of  0.25  inch  thickness  as  a  comparative  sample.  The  first  three 
entries  in  Table  I  list  values  of  R  along  with  the  maximum  and  minimum 
values  of  e'  for  the  SRRC  prototype  2.5  ns  air  line  system,  the  RG331  10  ns 
system  taken  prior  to  shipment  to  AFAL,  and  recen.,  results  taken  with  the 
AFAL  system.  As  expected,  the  percentage  ripple  of  the  2.5  ns  air  line  system 
is  much  smaller  than  the  AFAL  system. 

To  begin  a  meaningful  coaxial  line  investigation  before  the  arrival 
of  the  AFAL  equipment  at  SRRC  in  the  final  month  of  the  contract,  a  10  ns 
RG331  configuration  was  fabricated  at  SRRC.  The  configuration  was  first 
tested  with  a  2.5  ns  window  to  check  it  out  and  later  with  a  full  10  ns  win¬ 
dow.  Values  of  R  of  the  SRRC,  RG331  10  ns  configuration  for  both  the  2,5  and 

10  ns  window  are  given  in  Table  1. 

It  is  interesting  to  note  that  the  2.5  ns  RG331  system  displays 
less  ripple  than  the  10  ns  RG331  configuration,  although  the  delay  lines  are 
identical.  This  occurs  for  two  reasons.  First,  since  the  10  ns  window  gives 
four  times  as  many  frequency  points,  more  extreme  points  are  likely  to  occur. 
Secondly,  the  signal-to-noise  ratio  (SNR)  is  larger  for  the  noint  scan 
in  the  2.5  ns  window  than  the  1024  point  scan  in  the  10  ns  window.  The  other 

entries  in  the  table  will  be  referred  to  later. 

2.3  RIPPLE  SOURCE 

Tests  with  the  SRRC  10  ns  foam  cable  configuration  showed  the  same 
type  of  systematic  ripples  observed  in  the  AFAL  system.  An  examination  of  the 

’'•Registered  trademark  of  the  I.  E.  duPont  deNemours  and  Company. 
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TABLE  1 


COMPARISON  OF  MEASUREMENT  SYSTEM  RIPPLE 
ON  A  0.25  INCH  TEFLON  SAMPLE 

Time  Max  Min  Date  &  Place  of 


System 

Window 

2.5 

>  ns 

air 

line 

2.5 

ns 

10 

ns, 

331, 

AFAL 

10 

ns 

10 

ns. 

331, 

AFAL 

10 

ns 

10 

ns, 

331 

2.5 

ns 

10 

ns, 

331 

10 

ns 

10 

ns, 

air 

line 

10 

ns 

10 

ns, 

modified 

10 

ns 

f f  GHzl  fLGHzJ 


2.07 

2.0 

1.95 

9.6 

2.25 

9.2 

1.90 

8.1 

2.21 

9.6 

1.68 

6.1 

2.10 

9.6 

1.97 

2.0 

2.16 

9.2 

1.94 

8.3 

2.07 

0.9 

1.88 

8.6 

2.11 

2.0 

1.95 

6.7 

R 

Measurement 

6  % 

12/70 

SRRC 

17.5 % 

8/71 

SRRC 

26.5% 

6/72 

AFAL 

6.5% 

7/72 

SRRC 

11  % 

7/72 

SRRC 

9.5% 

7/72 

SRRC 

8  % 

9/72 

SRRC 
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delay  lines  using  TDR  methods  failed  to  reveal  any  discontinuities  in  the 
interior  of  the  RG331  cable.  However,  the  connectors  on  the  RG331  cable 
adjacent  to  the  sample  holder  did  exhibit  reflections  of  the  order  of  15  mV 
peak  to  peak  for  a  240  mV  ipput  step.  Further  experiments  with  one  of  these 
connectors  showed  a  reflection  coefficient  that  decreased  to  only  *21  dB  in 
the  vicinity  of  4  Gdz.  This  mismatch  appears  large  enough  to  be  a  major  con¬ 
tributor  to  the  observed  errors. 

As  another  check  to  determine  whether  the  connectors  were  the  main 
source  of  error  another  10  ns  configuration  was  constructed  .at  SRRC  using 
the  superior  GR900  air  line  series  connectors  with  two  5  ns  sections  of  air 
line.  This  system  furnished  two  10  ns  time  windows  but  an  extra  trans¬ 
mission  background  waveform  wi'h  the  sample  in  the  holder  had  to  be  included 
in  the  measurement  procedure  background  subtraction.  The  calculated  percen¬ 
tage  ripple,  R,  of  this  system  as  shown  in  Table  I  was  less  than  t.he  SRRC 
10  ns  foam  system  and  considerably  less  than  the  AFAL  10  ns  system.  These 
results  furnished  additional  evidence  that  the  connector  mismatch  was  the 
main  source  of  systematic  error.  r 


2.4  METHODS  OF  RIPPLE  REDUCTION 

From  the  evidence  of  the  previous  section,  it 'was  clear  that  any 
modification  which  resulted  in  a  lower  mismatch  ,  in  the  connectors  adjacent 
to  the  sample  holder  would  yield  reduced  systematic  ripples.  Three  modifica¬ 
tions  were  considered: 

a)  Replacing  the  existing  10  ns  and  5  ns  RG331  iines'by 
General  Radio  14  mm  air  lines  with  GR900  connectors. 

b)  Improving  the  design  of  the  RG331  to  GR900  connector  to 
reduce  the  mismatch, 

c)  Employing  a  software  correction  in  the  Fortran  program  for 
the  observed  connector  mismatch. 

The  first  is  an  obvious  solution  since  we  have  already  demonstrated 
that  an  air  line  system  produces  less  ripple.  .  However,  if  this  approach  were 
adopted,  the  10  foot  coiled  RG331  cable  inside  the  pulse  generator  box  would 
have  to  be  replaced  by  a  10  foot  run  of  air  line.  The  resultant  configuration 
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would  extend  for  15  feet  and  involve  an  undesirable  major  hardware  change. 

To  replace  the  cables  on  both  sides  of  the  sample  holder  by  air  line  also 
seems  like  aii  extreme  solution,  since  it  is  the  connectors  rather  than  the 
cables  themselves  that  cause  the  mismatch.  Therefore,  the  other  two  possi¬ 
bilities  were  carefully  examined  to  determine  whether  they  could  be  used  in 
place  of  the  complete  air  line  replacement. 

^  The  simplest  solution  appeared  to  be  to  improve  the  existing 

connector.  During  the  contract  to  build  the  prototype  material  measurement 
system,  the  RG331  to  GR900  connector  was  specially  designed  by  SRRC  since 
none  with  very  low  SWR  were  commercially  available.  The  original  and  modified 
transition  between  the  0.5”  foam  cable  and  0.56"  air  line  is  shown  in  Fig.  1. 
The  modifications  at  the  foam-air  interface  are  shov n  as  shaded  areas.  Al¬ 
though  the  region  throughout  the  transition  was  designed  to  be  50  n,  the  abrupt 
transition  at  the  foam-air  interface  contributes  an  equivalent  shunt  capacity 
which  produces  the  mismatch.  It  seems  however,  that  some  abrupt  transition 
must  always  be  present  between  the  two  geometries.  Therefore,  although  one 
must  accept  some  shunt  capacity,  the  match  can  be  improved  by  increasing  the 
series  inductance  at  the  discontinuity  to  compensate  for  the  existing  capa¬ 
city.  In  other  words,  since  the  shunt  capacity  reduces  the  characteristic 
impedance  ZQ  of  the  line  at  the  discontinuity,  the  geometry  should  be  changed 
to  increase  Zq  at  the  same  point.  In  modifying  the  connector,  ZQ  was  in¬ 
creased  by  reducing  the  step  at  the  interface  and  removing  some  of  the  foam 
dielectric  at  the  .junction.  The  compensation  was  done  in  steps  until  the 
amplitude  of  the  mismatch  as  observed  by  TDR  techniques  reached  a  minimum. 
Figures  2(a)  and  (b)  show  the  return  from  the  AFAL  end  line  connector  and  the 
modified  connector  on  the  10  ns  cable  adjacent  to  the  sample  holder,  res¬ 
pectively.  In  both  cases  the  input  was  a  250  mV  step.  Returns  examined  from 
other  unmodified  transitions  usually  showed  a  response  of  the  order  of  15  mV, 
The  photos  in  Fig.  2  demonstrate  the  improved  match  in  the  modified  connector. 
There  remains,  however,  a  residual  reflection  which  can  contribute  some 
systematic  error. 

The  third  approach  to  ripple  reduction  involves  measuring  the  mis¬ 
match  present  and  using  a  computer  correction  to  remove  the  effects  of  the 
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FIG.  1  Compensation  of  connector  shunt  capacitance  at  the  foam/air 
dielectric  interface. 
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reflection  from  the  final  result.  Therefore,  one  could  presumably  have  a 
poor  match  in  the  system,  as  long  as  it  was  known,  and  calculate  a  correction 
for  it.  A  discussion  of  this  problem  is  presented  in  Section  4,  since  the 
analysis  is  intimately  tied  in  with  calibration  with  standard  samples.  The 
advantage  of  this  method  is  that  even  if  the  line  mismatch  changes  over  a 
period  of  weeks  or  months,  the  procedure  will  measure  the  change  and  correct 
for  it.  Disadvantages  include  the  analysis  and  programing  effort  required 
to  achieve  an  operational  program,  the  additional  reference  waveform  required 
for  the  computer  correction,  and  the  possibility  of  inconclusive  results  for 
small  mismatches  because  of  contamination  by  random  noise. 

2.5  SELECTION,  FABRICATION  AND  TESTING  THE  MODIFICATIONS 

During  the  investigation,  each  of  the  three  above  methods  were 
studied  and  evaluated.  We  decided  that  a  combination  of  these  methods  was 
most  likely  to  produce  satisfactory  results.  First  the  connector  on  the 
10  ns  cable  adjacent  to  the  sample  holder  was  modified  to  compensate  for  the 
shunt  capacity.  This  approach  was  chosen  to  retain  the  original  pulse 
generator  configuration  with  the  coiled  cables  inside  the  box.  Secondly,  the 
5  ns  foam  end  line  was  replaced  by  an  air  line.  Unlike  air  line  replacement 
on  the  other  side  of  the  sample  holder,  this  did  not  involve  a  major  hard¬ 
ware  change  and  insured  an  excellent  match  on  the  end  line  side  of  the  holder. 
Special  care  was  taken  in  the  air  line  constructed  from  GR  precision  tube 
and  rod.  The  line  shown  in  Fig.  3  consists  of  2  long  and  1  short  section. 

The  middle  9  13/16"  section  has  GR900  AP  connectors  on  each  end.  Each  2654" 
section  has  GR900BT  connectors  on  each  end  which  provide  support  for  the 
inner  conductor.  The  upper  portion  of  the  crown  was  removed  from  one  of 
these  connectors  to  mate  with  the  sample  holder.  The  long  inner  conductors 
are  also  supported  at  the  center  of  each  line  by  a  J»"  polystyrene  disk.  To 
compensate  for  the  increase  in  capacitance  at  the  disk,  the  diameter  of  the 
center  conductor  was  reduced  as  illustrated  in  Fig.  3.  Experiments  showed 
negligible  reflections  from  the  completed  line.  The  air  line  is  supported 
by  an  I  beam  with  moveable  saddle  support  pieces.  Thirdly,  to  remove  the 
effects  of  any  residual  reflections,  a  computer  program  was  written  and  in¬ 
cluded  in  a  modified  Fortran  routine  for  the  measurement  sequence. 


A 


CROWN  REMOVED 


POLYSTYRENE  FOAM 
SUPPORT 

I  mmm 


.CONNECTOR  DIAMETER  UNDERCUT 
0.009" 

tzzzn 


viirnriji 


H  }*~  1/4" 


g  TYPE  900 BT  CONNECTORS 
□  TYPE  900  AP  CONNECTORS 


CENTER  CONDUCTOR  JOINT  AND 
SUPPORT  IN  THE  MIDDLE  OF  THE 
26 1/4"  SECTIONS. 


FIG.  3  GR  900  air  line  configuration. 


The  above  modifications  were  made  on  part  of  the  AFAL  system 
returned  to  SRRC  in  mid-August  and  resulted  in  a  reduction  of  systematic 
error.  Since  the  connector  modification  reduced  the  mismatch  to  a  very  low 
level,  no  significant  improvement  was  noted  when  the  computer  correction  was 
applied  and,  therefore,  the  data  in  the  remainder  of  this  section  is  shown 
without  software  correction.  Figure  4  shows  a  comparison  of  the  complex  p, 
and  e  for  the  original  system  (Fig.  55,  Ref.  1)  and  the  improved  system. 

Note  the  absence  of  the  large  variations  around  3  GHz  originally  present. 

The  poorer  low  frequency  values  for  the  improved  system  are  random  variations 
due  to  noise.  The  percentage  ripple  R  computed  for  the  modified  system  (which  on 
different  days  varied  between  7.5%  and  9%  -  the  most  common  value  of  8%)  is 
included  in  Table  1.  Another  before  and  after  comparison  is  illustrated  in 
Fig.  5  where  the  real  part  of  e  for  nylen  and  a  ferrite  is  plotted  for  the 
AFAL  system  after  it  was  returned  to  SRRC  and  the  improved  system.  The  AFAL 
equipment  returned  to  SRRC  had  an  obvious  discontinuity  in  the  outer  con¬ 
ductor  of  a  connector  near  the  sample  holder.  The  data  in  Fig.  4  was  taken 
after  the  connector  was  realigned  and  mating  surfaces  cleaned.  These  com¬ 
parisons  between  the  improved  and  original  system  illustrate  the  reduction  in 
systematic  error  achieved  in  the  coaxial  line  investigation. 

Several  improvements,  not  required  by  the  contract,  were  made  on  the 
components  in  the  pulse  generator  box  of  the  AFAL  system.  The  avalanche 
transistor  module  was  rebuilt  to  have  superior  physical  design.  The  lantern 
battery  bias  supply  for  the  step  recovery  diodes  was  replaced  by  a  dc  power 
supply,  and  a  fuse  was  installed  in  the  power  line.  Figure  6  now  shows  the 
new  power  supply  block  diagram.  The  RG331  foam  line  between  the  sample  head 
and  sample  holder  was  replaced  by  another  RG331  line  whose  impedance  was 
closer  to  50  ohms. 
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FIG.  5  Comparison  of  e'  results  for  original  and  modified  system 
for  nylon  and  ferrite  JF-1, 
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SECTION  3 

LOW  FREQUENCY  RESPONSE 


3.1  INTRODUCTION 

Dsta  below  about  500  MHz  usually  appears  more  erratic  than  results 
at  high  frequencies.  Some  reasons  for  these  randomly  occurring  values  could 
be  in  setting  the  reference  points,  positioning  the  sample, or  noise.  Since 
it  has  been  noted  that  repetitive  runs  on  the  same  undisturbed  sample  also 
display  erratic  low  frequency  data  points,  the  most  probable  cause  is  random 
noise. 


The  effects  of  random  noise  in  and  Sgj  are  more  serious  at  the 
lowest  frequencies  where  is  near  0  and  is  near  1.  The  reason  for 

this  is  that  since  the  calculations  of  p,  and  s  depend  essentially  on 
and  1-S0j,  the  same  additive  errors  in  these  quantities  produce  much  larger 
percentage  errors  when  S^.  is  near  0  and  is  near  1.  This  can  be 
illustrated  by  a  simple  example  using  experimental  values  of  lsjjl  and  I S2j l • 

TABLE  2 

SCATTERING  PARAMETER  ERRORS 


FP.EQ. 

V 

tS2l  1 

ASSUMED 

ADDITIVE 

ERROR 

%  ERROR 

IV 

%  ERROR 

i“lS2li 

100  MHz 

0.009 

0.995 

0.001 

n% 

2Q$ 

3000  MHz 

0.196 

0.974 

0.001 

0,5$ 

4% 

Thus  oue  would  certainly  expect  more  error  at  100  than  3000  MHz.  Of  course 
a  more  complete  error  analysis  would  have  to  account  for  phase  errors  also, 
but  the  above  illustrates  the  main  point. 

3.2  EXPERIMENTS  AND  RESULTS 

The  previous  example  suggests  that  the  low  frequency  behavior  can 
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be  improved  if  the  SNR  is  increased.  However,  the  price  paid  for  better  SNR 
at  lower  frequencies  is  usually  poorer  SNR  at  higher  frequencies.  One  method 
to  increase  the  low  frequency  SNR  is  to  use  a  material  with  a  large  dielectric 
constant  as  a  test  sample  to  cheek  the  system.  One  notes,  for  example,  that 
ferrite  samples  usually  exhibit  a  fairly  flat  low  frequency  response.  This 
method  is,  of  course,  inapplicable  if  the  low  frequency  behavior  of  a 
particular  material,  such  as  Teflon  is  desired.  A  wider  excitation  pulse 
will  also  furnish  a  hight-r  low  frequency  SNR,  An  experiment  was  performed 
with  a  250  ps  pulse  which  had  about.  12  dB  more  spectral  amplitude 
than  the  60  ps  pulse  at  low  frequencies.  As  expected,  the  data  showed  a 
smoother  lower  frequency  response  with  much  more  erratic  high  frequency 
values.  To  use  this  method  would  require  another  pulse  generator  or  wave¬ 
form  integrator  in  the  system.  Perhaps  the  simplest  approach  to  improve  low 
frequency  behavior  is  the  use  of  a  thicker  sample  that  will  give  more 
reflection  and  less  transmission  than  a  shorter  sample.  However,  in  addition 
to  the  poorer  high  frequency  SNR  a  longer  sample  can  produce  a  singularity  in 
p,  and  e  when  the  sample  thickness  approaches  a  half  wavelength.  It  has 
also  been  observed  that  after  the  singularity,  the  values  of  u,  and  e  are 
grossly  incorrect.  Since  the  longer  sample  approach  appeared  most  advanta¬ 
geous,  more  analysis  was  done  to  see  if  the  response  around  and  after  the 
discontinuity  could  be  improved. 

3.3  THE  THICK  SAMPLE  PROBLEM 

The  values  of  p,  and  e  are  eri  stic  in  the  neighborhood  of  the 
singularity  for  the  same  reason  that  they  a  t-  erratic  at  low  frequencies. 

Namely,  S  ^  is  near  0  and  S2^  is  close  to  1.  In  fact,  the  response  is 
usually  worse  since  the  SNR  of  the  incident  signals  is  usually  lower  near 
the  higher  frequency  singularity  than  the  one  at  zero  frequency.  The  singularity 
cannot  be  avoided  with  the  present  program  since  division  by  S^  s»  0  occurs. 
Although  'iata  in  the  region  around  the  singularity  will  continue  to  be 
erratic,  one  should  bo  able  to  improve  the  values  of  p,  and  e  above  the 
singular  frequency. 

An  examination  of  the  present  program  revealed  that  p,  and  e  above  the 
singularity  were  incorrect  because  the  argument  of  the  quantity  z  =  exp 
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C-ju/lA  '«  'd/c),  where  d  is  the  sample  thickness,  (see  also  Eq.  (2),  p,  3, 

Ref,  2),  hereafter  referred  to  as  arg  z,  was  not  computed  correctly  for  arg  z 
less  than  -it  .  If  arg  z  is  wrong,  Cg  (Eq.  (9),  p.5,  Ref.  2)  and  therefore 
p/c'  will  also  be  incorrect.  A  study  of  Eqs.  (2)  through  (6)  in  Ref.  2  also 
showed  that  the  singularity  occurs  when  arg  z  =  -n.  The  incorrect  values 
occurred  because  the  arctangent  routine  that  computed  arg  z  only  yielded  angles 
between  ±n.  Thus  for  z  in  the  third  quadrant,  the  program  gave  arg  z  -  a, 
whereas  the  correct  ary  z  should  have  been  a  -  .  To  correct  this,  a 

change  was  made  in  the  subroutine  CPLEX  to  give  urg  z  between  0  and  -2n.  The 
modification  is  listed  in  Appendix  B.  This  revised  program  was  tried  and 
yielded  the  correct  values  of  arg  z.  Figure  7(a)  and  (b)  offers  the  com¬ 
parison  between  s'  for  a  0.499  inch  Teflon  sample  computed  with  the  original 
program  and  the  modified  program.  Note  that  the  region  of  grossly  incorrect 
values  for  t'  extends  for  a  considerable  frequency  interval  on  both  sides  of 
the  singularity  at  8.4  GHz.  This  is  probably  due  to  the  decreased  SNR  of 
l-S^j  for  the  longer  sample  at  higher  frequencies.  It  should  be  mentioned 
that  the  revised  subroutine  will  still  give  incorrect  angles  for  arg  z  <  -  2tt. 

A  more  involved  program  change  is  required  to  track  the  phase  beyond  one  com¬ 
plete  revolution. 

If  very  accurate  data  is  desired  at  both  low  and  high  frequencies, 
one  could  test  both  a  thin  and  thick  sample  of  the  same  material.  The  low 
frequency  data  would  be  furnished  by  the  thick  sample,  while  the  thinner 
sample  would  supply  the  higher  frequency  results. 
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FIG.  7  Angle  ambiguity  of  e#  for  a  0.499  inch  Teflon  sample 
computed  with 

a)  Original  program 

b)  Modified  program 
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SECTION  4  ,  , 

CALIBRATION  WITH  STANDARD  SAMPLES 

4.1  INTRODUCTION 

Another  approach  investigated  to  reduce  the  systematic  ripples  was 
based  on  software  correction.  For  example,  'it  should  be  possible  to  employ 
a  standard  sample  and  deduce  a  system  correction  factor  by  Comparing  the 

i 

measured  to  known  values.  Although  computer  techniques  appeared  feasible, 
our  plan  of  attack  was  to  first  reduce  mismatches  as  much  as  possible  via 
hardware  improvements,  i.e,  connector  modification  and  line  replacement, 
and  then  to  try  the  software  correction  on  any  remaining  residual  reflections^ 
This  plan  was  adopted  because  we  felt  in  the  short  contract  period,  the  hard¬ 
ware  modifications  showed  the  greatest  probability  of  success.  As  it  .turned  out 
after  the  hardware  modifications  reduced,  the  ripple  to  a  lower  value,  the 

1  i 

subsequent  software  correction  did  not  significantly  improve  the  results. 

4.2  ANALYSIS  ,  . 


In  order  to  calculate  a:  correction  factor,  one  has  to  determine  the 
relationship  between  the  measured  and  true  values.  One  of  the  -simplest  models 
for  line  mismatch  is  shown  in  Fig.  8(a)  for  a  single  shunt  discontinuity  with 
reflection  coefficient  r(<u)  a  distance  d  to  the  left  of  the  sample. 
Fortunately  this  model  should  be  a  good  approximation  to  the  modified  system  • 
where  the  air  line  to  the  right  of  the  sample  contributes  almost 'no  discon¬ 
tinuity  while  the  modified  connector  to  the  left  of  the  sample  contributes  a 

small  shunt  discontinuity.  The  bounce  diagram  generated  by  this  configuration 

b  1  ! 

is  shown  in  Fig.  8(b),  where  VR(u>)  and  V^.(cu)  are  the  reflected  and  transmitted 
signals  respectively  at  a  particular  frequency,  and,  S^(u>)  and  ate 

the  true  values  of  the  scattering  parameters.  Summing’all  the  reflected 


FIG.  8  Reflection  and  transmission  in  the  presence  of  a 
discontinuity. 


where  z  =  exp  (-2j<ud/v)  and  vis  the  propagation  velocity  in  the  line.  The 
first  reflection,  r,  is  not  included  since  it  arises  from  a  part  of  the  time 
waveform  before  the  reflection  time  window.  Similarly  the  sum  of  the  trans¬ 
mitted  phasors  are 
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In  the  experimental  procedure,  the  measured  value  of  the  scattering 
m 

coefficient,  S  ,  is  found  from  ratio 
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where  the  background  signal,  VRR,  is  found  from  VR  with  =  0,  giving 

VRB  =  0 


(4) 


The  incident  waveform,  VRI,  is  found  by  placing  a  short  at  the  sample  position, 
or  equivalently  setting  =  -1.  The  signal  measured  is  the  negative  of  the 
incident  phasor.  Therefore, 


v  =  il  + 

RI  1  +  rz 


(5) 


Substituting  (1),  (4)  and  (5)  into  (3)  gives 


-  c  l+o 

i  r  5n  i  -  ps 


(6) 


11 


where  we  have  defined  p  =  rz  and  call  it  "the  correction  factor".  Actually 
p  is  simply  the  reflection  coefficient  looking  into  the  line  at  the  sample 
position. 

Similarly  the  measured  value  of  the  scattering  parameter,  S^j,  is 
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found  experimentally  from 


,m  _  VT  ~  VTB 
2!  vTj  -  VTB 


(7) 


The  background  wave  in  the  transmission  window  is  found  by  removing  the 
sample  and  terminating  the  line,  i.e.  setting  =  0,  giving 

VTI  =  0  .  (8) 

The  incident  signal,  V^,j  is  measured  with  the  sample  removed  or  equivalently 
setting  S2j=l  and  Sjj=0  yielding 

VTI  =  (1+r)  (9) 


Substituting  Eqs.  (2),  (8),  and  (9)  into  (7)  gives 

cm  -  S21 

21  l-pSu 


(10) 


The  result  is  that  Eqs.  (6)  and  (10)  provide  the  desired  relations  between 
the  measured  and  actual  values  of  the  scattering  coefficients  in  terms  of  the 
correction  factor  p  .  It  should  be  mentioned  that  although  the  expressions 
were  derived  for  a  single  shunt  discontinuity  they  are  also  valid  for  any 
line  reflections  describable  by  a  reflection  scattering  parameter  A^f^), 
simply  by  replacing  p  by  A^.  The  correction  factor  p  can  be  found 
from  known  and  measured  values  of  a  standard  sample  by  rewriting  Eq.  (6)  as 


P  - 


S'"  -  S 
-  11  11 


11 


(l+s^; 


(id 


After  p  has  been  determined  from  the  standard  sample  for  each  frequency,  it 
can  then  be  used  to  correct  the  measured  values  and  of  any  sample. 


4.3  MEASUREMENTS  USING  A  STANDARD  SAMPLE 


A  computer  program  was  written  so  that  the  correction  factor  p 
was  computed  from  a  standard  sample  using  Eq.  (11)  and  subsequent  corrections 
applied  to  other  samples.  A  0.25  inch  thick  Teflon  sample  was  used  as  a 
standard  with  an  assumed  e  =  2.03  +  jO  and  p,  =  1  +  jO.  The  actual  values 
of  S^  and  S£j  for  Teflon  were  computed  from  Eqs.  (1)  through  (4)  of  Ref. 

2.  The  results  were  somewhat  disappointing  since  the  corrected  data  did  not 
display  significantly  less  ripple  than  the  uncorr&eted  data.  Several  reasons 
can  be  advanced  for  this  behavior.  In  the  first  place,  the  experiments  were 
performed  with  configurations  exhibiting  quite  low  reflections.  Thus  there 
was  little  error  to  correct  and  any  corrections  that  did  survive  could  suffer 
serious  noise  contamination.  Secondly  we  had  to  assume  an  e  and  p,  for  the 
standard  sample.  If  the  choice  was  incorrect,  it  might  contribute  more  error 
than  was  originally  present.  Because  of  these  problems  another  method  was 
sought  which  might  provide  better  results. 

4.4  THE  DISPLACED  SHORT  METHOD 

In  the  previous  analysis  a  standard  sample  is  used  to  find  a 
correction  factor  for  the  line.  However,  to  characterize  the  line,  a  stan¬ 
dard  sample  is  not  necessary  since  two  measurements  on  a 'short  and  displaced 
short  will  accomplish  the  same  task.  The  displaced  short  also  minimizes 
several  disadvantages  of  the  previous  approach.  In  the  first  place  the  SNR 
of  the  reflected  wave  is  higher  since  the  signal  is  larger.  Secondly,  no 
assumptions  are  necessary  about  the  short's  parameters  except  that  its 
reflection  coefficient  is  -1. 


The  relation  between  p  and  displaced  short  measurements  can  be 
found  from  the  previous  equations.  If  the  short  is  a  distance  dj  from  the 
discontinuity,  the  reflected  wave  is  given  by  an  equation  similar  to  (5), 


(l+r)^z 

V  _ _ i. 

R1  l+rzj 


(12) 


where  Zj  =  exp(-2j<udj/v).  If  the  short  is  now  displaced  an  additional 
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distance  D  to  the  right,  so  its  new  distance,  d2  from  the  discontinuity 
is  given  by 


d2  =  dj-  +  D 

the  reflected  wave  is  now 

(l+r)2z9 

i;  =  _  - i 

R2  l+rz2 


The  ratio  of  the  two  reflections,  V,  is 


V  z.U+rzJ  z2 


(13) 


(14) 


where  as  before,  p  =  rZj.  Noting  Zj/z2  =  exp(2ji»D/v)  the  correction  factor 
can  be  found  from  the  measured  ratio  and  the  calculated  value  of  Zj/z.,  as 


2iupD 

V-e 

P  "  “  V-l 


(15) 


Another  program  was  written  to  compute  p  from  Eq.  (15)  rather 
than  (11).  To  obtain  the  necessary  data  on  the  displaced  short,  VR2  , 
another  reference  waveform  was  acquired.  Experiments  using  the  displaced 
short  method,  as  in  the  calibrated  sample  method,  did  not  show  significant 
variations  between  the  corrected  and  uncorrected  values.  Figure  9  illustrates 
values  obtained  for  e'  with  and  without  the  displaced  short  correction  for 
nylon  and  ferrite.  To  study  possible  differences,  comparisons  were  made  on 
data  taken  on  different  day*  to  separate  random  from  deterministic  errors. 

A  careful  examination  revealed  that  some  of  the  systematic  error  was  reduced 
by  the  correction,  while  other  apparently  systematic  errors  remained.  We 
hypothesized  that  the  residual  mismatch  is  so  small  that  software  correction 
does  not  make  much  difference. 
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FIG.  9  Comparison  of  data  with  and  without  displaced  short  corrections. 


-  26  - 


To  determine  whether  the  correction  was  working  at  all,  we  purposely 
introduced  a  large  mismatch  in  the  line  by  inserting  a  Teflon  disk  against 
the  sample  holder  connector.  The  results  showed  the  displaced  short  pro¬ 
cedure  was  indeed  reducing  the  systematic  error.  For  example,  a  measurement 
on  nylon  with  the  uncorrected  data  gave  a  percentage  ripple,  R  =  33%,  whereas 
the  corrected  data  gave  R  =  17%.  It  should  be  stressed  that  efforts  to 
reduce  the  systematic  error  do  not  reduce  the  random  error.  Therefore,  the 
erratic  results  at  lower  frequencies,  which  are  due  to  noise,  are  not.  improved 
by  either  the  calibrated  sample  or  displaced  short  procedure. 

4.5  CONCLUSIONS  AND  OPTIONS 

In  conclusion,  we  have  demonstrated  the  feasibility  of  reducing 
systematic  error  by  computer  correction.  However,  since  the  mismatch  was 
already  reduced  to  a  low  level  by  hardware  corrections,  additional  improvement 
by  software  techniques  in  the  modified  system  was  nearly  imperceptible.  Of 
the  two  correction  methods  developed,  the  displaced  short  seems  superior 
because  of  the  larger  SNR  of  its  waveforms  and  the  fact  that  no  constitutive 
parameters  are  required. 

In  Appendix  A,  we  have  included  a  revised  main  program  to  compute 
n  and  e  with  the  displaced  short  procedure.  The  operator  has  the  option 
of  scanning  the  waveforms  as  described  in  Ref.  1  and  using  the  original 
Fortran  program  or  including  the  displaced  short  waveform,  with  the  revised 
program.  The  only  change  in  the  measurement  procedure  is  the  acquisition  of 
another  reference  waveform  after  the  incident  waveform  in  the  R  window  is 
scanned.  The  revised  procedure  below  should  be  inserted  after  the  B6  operat¬ 
ing  sequence  as  given  on  p.  18  of  Ref.  1  and  before  step  B7. 

Displaced  Short  Waveform  Modification 
B6A  COAXCR)  is  lit 

The  metal  short  already  in  place  in  the  coaxial  line  1  inch 
from  the  incident  connector  should  now  be  positioned  with  the 
supplied  gauge  block  1.4  inches  from  the  incident  connector. 
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To  continue,  push  SHORT 

Data  from  the  scan  in  step  B6  will  first  be  written  on  tape  if 
WRITE  (red)  is  lit,  and  SHORT  REFL  (red)  will  again  be  lit 
during  the  scan. 

B7  COAX(R)  is  lit 

Continue  original  operating  sequence 

We  have  demonstrated  that  software  corrections  are  feasible,  but  do 
not  significantly  improve  results  if  the  residual  mismatch  is  small. 
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SECTION  5 

LOW  LOSS  DIELECTRIC  MEASUREMENTS 

5.1  INTRODUCTION 

A  secondary  task  during  this  contract  was  the  investigation  of 

techniques  for  adapting  the  AFAL  measurement  system  to  low  loss  dielectric 

-3 

materials.  Materials  with  loss  tangent  as  sm3ll  as  10  are  of  interest, 
particularly  over  the  frequency  range  1  to  10  GHz.  The  prototype  AFAL  system 
was  designed  to  measure  e  and  p  with  errors  of  the  order  of  *0,1,  Al¬ 
though  these  errors  are  acceptable  for  testing  high  loss  RAM  type  materials, 
they  are  too  large  for  specifying  losses  in  radome  type  materials.  Low  loss 
measurements  are  inaccurate  because  the  loss  introduced  by  the  nominal  0.25 
inch  samples  is  smaller  than  the  measurement  uncertainty.  For  example,  a 
signal  transmitted  through  a  0.25  inch  sample  of  a  material  with  6*  =  2  and 
tan  6  =  10“^  is  attenuated  by  about  1  part  in  104  at  1  GHz.  If  the  incident 

signal  has  an  SNR  of  60  dB,  which  translates  into  a  measurement  uncertainty 

3 

of  approximately  1  part  in  10  ,  then  the  uncertainty  is  10  times  larger  than 
the  quantity  we  seek  to  determine.  Obviously  no  meaningful  measurement  can 
be  made  in  this  situation. 

5.2  THE  LONG  SAMPLE  METHOD 

One  obvious  way  to  increase  measurement  accuracy  is  to  introduce 
more  transmission  loss  by  using  a  longer  sample.  A  convenient  geometry  for 
time  domain  measurements  is  a  sample  of  material  completely  filling  the  trans¬ 
mission  line  backed  by  a  short  circuit.  For  this  configuration  the  first 
reflection  from  the  air-dielect.ric  interface  is  used  to  determine  6 '  , 

while  the  second  reflection,  which  traverses  the  sample  twice,  gives  e"  . 

The  relation  between  the  sample  attenuation  and  measurement 
uncertainty  can  be  found  from  the  following  expressions.  The  total  dielectric 
loss,  or  ,  for  the  TEM  propagation  mode  is  given  in  dB  by  Ref.  (3), 

[a]  =  4.62  FdA'  tan  6  (16) 

where  F  is  in  GHz,  d  is  the  sample  length  in  inches  and  the  notation  is 
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introduced  where  square  brackets  around  a  quantity  denote  dB,  i.e. 

[a]  =  20  log^  a.  The  transmission  loss  of  the  material  is  computed  from  the 
ratio  of  the  Fourier  transform  of  the  transmitted  incident  waveform.  These 
waveforms  are  contaminated  with  noise.  If  gaussian  noise  is  assured,  the 
standard  deviation,  a  ,  of  the  error  is  given  in  dB  by  Ref.  (4) 

[a]  a  2*1  a  +  a2)^  (17) 

bi 

where  S.  is  the  SNR  of  the  incident  wave.  Note  that  both  and  a  are 
expressed  in  Eq,  (17)  as  voltage  ratios.  Equation  (17)  is  a  valid  approxi¬ 
mation  for  S./a  >  10.  The  criterion  for  accurate  measurements  is 
[or]  »[o]. 

Equations  (16)  and  (17)  can  be  used  along  with  the  incident  SNR  to 

determine  what  sample  lengths  are  required  for  good  measurements.  For 

example,  suppose  one  desires  to  measure  a  dielectric  with  s'  =  2  and  tan  6  = 
-3 

10  .  What  minimum  length  is  required  to  measure  this  material  over  the 
frequency  range  1  to  10  GHz,  such  that  [a]  >  10  [c]  over  the  entire  range? 

The  problem  can  be  solved  graphically  using  the  experimentally  determined 
incident  SNR  found  in  Table  3. 


TABLE  3 

INCIDENT  SNR  FOR  VARIOUS  FREQUENCIES 

Freq.  (GHz)  1  3  5  8  10 

Incident  SNR  (dB)  62  59  54  48  37 

-3 

Since  a  will  be  close  to  1  for  a  tan  6  =  10  material,  Eq.  (17)  can  be 
approximated  by  [a]  as  12.3/S.  .  This  quantity  is  plotted  in  Fig.  10  along 
with  10  [a].  As  a  first  trial,  a  length  d  =  12"  is  chosen,  and  [a]  com¬ 
puted  from  Eq.  (16)  and  plotted.  The  graph  shows  that  the  inequality 
[a]  >  10  [a]  is  not  satisfied  at  the  low  and  high  frequency  end.  A  new  d 
of  27"  is  determined  by  setting  [a]  equal  to  10  [a]  at  10  GHz.  As  shown  in 
Fig,  10,  the  27"  length  satisfies  the  inequality  over  the  entire  frequency 
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FIG,  10  Graph  for  determination  of  minimum  sample  length. 
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range  and  therefore  is  the  minimum  length  sought. 

Thus  far  the  calculations  have  not  accounted  for  copper  loss  in  the 
portion  of  the  line  occupied  by  the  sample.  Since  this  loss  depends  on  the 
dielectric  constant  in  the  line,  it  does  not  cancel  out.  When  the  copper 
loss  approaches  the  order  of  the  dielectric  loss,  it  must  be  accounted  for. 
Perhaps  the  most  practical  method  of  achieving  this  is  through  a  computer 
correction  of  the  measured  data. 

A  variation  of  the  long  sample  method,  based  on  measurements  of  the 
reflections  making  four  or  six  traverses  of  the  sample,  should  be  considered 
if  the  required  sample  length  for  the  two  traverse  reflections  becomes  un¬ 
wieldy.  Sources  of  difficulty  for  this  multiple  pass  method  will  be  in 
making  measurements  where  the  incident  pulses  may  not  be  in  the  same  time 
window  as  the  transmitted  pulse  and  accounting  for  copper  losses. 

Strip  line  sample  holders  do  not  seem  to  offer  any  advantages  over 
the  coaxial  line  unless  the  material  is  only  available  in  thin  sheet  stock. 

Even  then  it  might  be  preferable  to  stamp  disks  from  the  sheet  and  mount  a 
row  of  disks  in  the  coaxial  line.  The  higher  order  modes  introduced  at  the 
stripline  to  coaxial  junction  will  complicate  the  measurement  interpretation. 

In  conclusion,  the  long  sample  method  appears  to  be  the  best  way 
at  present  to  measure  low  loss  materials.  A  long  sample  holder  backed  by  a 
short  circuit  could  be  made  to  accommodate  the  sample.  The  measurement  system 

should  be  versatile  enough  so  that  different  sample  lengths  can  be  used. 

—  •> 

That  is,  longer  samples  will  be  required  for  materials  with  tan  6  =  10~‘'  than 
for  tan  6  =  10-1.  The  equations  given  in  this  section  can  be  used  to  deter¬ 
mine  the  lengths  required.  This  versatility  can  be  included  by  having  variable 
position  reference  waveforms  and  variable  time  windows.  It  is  suggested 
that  another  network  configuration,  compatible  with  the  existing  instrument 
rack,  be  built  to  perform  the  low  loss  measurements. 
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APPENDIX  A 


REVISED  MAIN  PROGRAM  WITH  DISPLACED  SHORT  CORRECTION  ' 

i 

*  »  ! 

The  new  main  program  is  listed  below.  The  data  card  with  the  dis 

placed  short  distance  of  400  mils  should  be  inserted  immediately  before  the 

i 

sample  thickness  cards. 
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CALL  CRLEXtUPR.URI.DNR.il. RHOR  .RHOI  .21 

c 

REF 

0010TZ 

CALL  CRIEXtTC  1 1. 11  .TCI  Z.I1  .TAt  1  .Il.TA  tZ. I1.REFR.REFI.  Zl 

c 

CHINGE  SIN  REF 

001 1 27 

REFR3-REFR 

001135 

REFIr-REFl 

c 

CORRECT  NITH  RHO 

0011IZ 

L ALL  CPLEXtt ..REFR.REFI.RMOR  *RMOI  .AR.AI.il 

001 1S1 

CALL  CRL EXtRSFR .REF I. 1 ..AR.AI.REFR. REF I.Z1 

001Z01 

RMOD  3SIRTfREFR.REFR.REFI.REFI 1 

ooizzs 

RPH3ATANZtREFI.REFRl.FC 

c 

TRA 

001 Z3T 

CALL  CRL EXtTDtl.il .TDtZ.Il.TBtl.il, TBtZ.Il.BR. 81. 21 

001ZT3 

CALL  CRLEXlBR.8IfDR.DI.AR, AT, 3| 

001 301 

RHI3FACA.F 
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lot 

001313 

CALL  CPI  EX!  AR  •  Al.COSt  PHI1  SIN  (PHI  t.TRAR  .THAI 

JOS 

00133T 

CALL  CPLEX fREFR.REF I <RHOR  .RHOI  .AR.AI.lt 

IBS 

001350 

AR=t.-AR 

LOT 

OOtlST 

AI=-AI 

IBS 

OOI3S5 

CALL  CPlEXfTRAR.TRAI.An.AI.TRAR.TRAI.t  ) 

IOJ 

00 1  STS 

TROD :SSRT(TRAR*TRAR*TRAI«TRAII 

no 

001 1 22 

TPH=ATAH2tTRAI.TRARl»FC 

tit 

C 

TI.S2 

112 

oottss 

ftR-TRAR.REFR 

113 

OOlttl 

»ii=trai«ref: 

tit 

OOlttT 

ar=trar-refr 

IIS 

oottss 

ai=trai-refi 

IIS 

C 

X 

I  IT 

oottss 

CALL  CPLEXf9tR.flT.ARtAr.DR.Dr.lt 

IIS 

OOltTt 

AR=1.-0R 

IIS 

001503 

AI--DI 

120 

001510 

CALL  CPI  EXfAR.AX.2.*REFR<  2  .  .REF  TtXR.XT.2l 

121 

C 

GANNA 

122 

001535 

CALL  CPLEXtXR.XI.XR.XT.BR.BI.lt 

123 

ootsts 

BR=BR-1 . 

I  2t 

ootsss 

CALL  CPLEXfBR.BT.DR.0T.AR.AI.3t 

125 

ootsss 

8R-XR.AR 

12S 

001  STS 

Gt=XT«AI 

I2T 

001S03 

IF  fSBRT fGR.GR.BT.GI 1  .LE.  1 « t  GO  To  9 

121 

ootsss 

GR-XR-AR 

129 

ootstt 

GT-XI-AI 

130 

001 652 

9  CONTINUE 

l'Jl 

C 

2 

132 

133 

ootsss 

CALL  CPLEXtftR.VlI.GR. GI.AR.AI.lt 

CALL  CPLEXI91R-GR.fi I-GI.l.-AR.-AI»2R«2I.2t 

I3t 

C 

Cl 

13S 

001T25 

CALL  CPLEX  f  |  . «CR  >GI  .1.-0R.-GI.C1R.CII.2I 

t3S 

OOITST 

CALL  CPLEX f2R.2I.DUN.DUH.AR. AT, t| 

13T 

BOITTO 

C2RSAR*FACB/F 

t  31 

002001 

C2I=AI.FAC8/F 

139 

002011 

CALL  CPLEXtClR.ClI.-C2I.C2R.XHR.XHT. I  1 

ItO 

00202T 

CALL  CPLEXf-C2t.C2R.ClR. Cl I.EPR.EPI.Xt 

ttl 

0G20SS 

HRTTEfS.SOtF.XHR.XHl.EPR.epI 

1 12 

0020S2 

..RHOO.RPH.THOO.TPH 

MS 

002052 

..2R.2X.GR.GI 

Itt 

002! to 

GO  FORHATfF  10. 0.  1 2F 1 0. tt 

its 

002151 

XfNtlF 

It’S 

OOZISt 

TfNtxEPR 

1ST 

002ITT 

F=F.OELF 

1 1 T 

•DIAGNOSTIC* 

THIS  statehert  HUHBER  IS  DOUBLY  oefined 

its 

00220S 

11  CONTINUE 

1 1 9 

00221 t 

Nf ALIN2-N 1 ♦ 1 

150 

002221 

CAU  PLOTfl.XfNtl.YfNtt.NfAL.Jt 

1st 

0022t2 

J=J*I 

152 

0022tS 

READ  15.101  HORE 

1  S3 

0022S0 

IFIHORE.GT.OI  GO  TO  10! 

t$t 

002270 

toz  CONTINUE 

tss 

0022T0 

STOP 

ISS 

0022T3 

tot  WRITE  fS .tOOt 

1ST 

. 002301 

too  FORHATf*  FATAL  ERROR*! 

ISS 

002312 

GO  TC  t02 

159 

00231 1 

end 
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APPENDIX  B 


REVISED  CPLEX  SUBROUTINE 

In  order  to  remove  the  incorrect  angle  computation  as  described  in 
Section  4,  the  end  of  the  CPLEX  Subroutine  in  the  original  program  should  be 
modified  as  given  below. 

End  of  Original  CPLEX  Subroutine 
5  AR  =  ALOG(R) 

AI  =  TH 
GO  to  2 
END 

Replace  above  statements  with 

5  AR  =  ALOG(R) 

AI  =  TH 

IF  (AI.  GT.  O.O)  GO  TO  6 
GO  to  2 

6  AI  =  AI  -  6.28318 
GO  to  2 

END 


